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Abstract. The structure of the ground state beyond the instability point of the 
Q> . quasiparticle system with Fermi-step momentum distribution is studied within 

ON I the model of a Fermi liquid with a strong repulsive interaction. A ground state 

rearrangement occurs as the interaction strength is increased beyond a defi- 
^ ■ nite critical value. Numerical investigation of the initial stage of this structural 

transition shows that there are two temperature regions, corresponding to differ- 
| ent scenarios of the rearrangement. While for temperature T larger than some 

characteristic temperature Tq the behaviour of the system is the same as that in 
t— I , the case of the fermion condensation, for T<Tq the intermediate structure with 

multiconnected quasiparticle momentum distribution arises. The transition of 

■ this structure to the fermion condensate at increasing interaction strength is 
. discussed. 

■ 

o ■ 

■ 1. Introduction 

The question of an applicability of the Landau Fermi liquid theory |I| for describing 
properties of strongly correlated Fermi systems has been discussed since a long time. It 
is known that this theory is not valid for one dimensional (ID) systems 0. For such 
systems, the concept of Luttinger liquid || with a single-particle Green function con- 
taining no quasiparticle pole is usually introduced, instead of the quasiparticle picture. 
The frontiers of the non- Fermi liquid view cover also strongly correlated 2D liquids j|, [J, 
since HTC materials with quasi-2D structure possesse properties in contradiction with the 
expected ones in the Landau theory. However the recently measured electronic spectra 
of such materials [§, [7|, || show evidence for the presence of a quasiparticle pole in the 
single-electron propagator. At the same time, new possibilities were found within the 
quasiparticle approach in refs. [9-12], where quasiparticles with momentum distributions 
differing from the ones assumed in the Landau theory are introduced. This new class of 



systems with the presence of a fermion condensate, as predicted in refs. [|TT| |T2|| , possesses 
a rich variety of properties [12-15], some of which are characteristics of a non- Fermi liq- 
uid behaviour. As discussed within different models [12-14, 16], a state with the fermion 
condensate arises as a result of the rearrangement of the ground state of the system. 
This rearrangement in a system of quasiparticles, whose momentum distribution has the 
shape of a Fermi sphere with an occupation number slightly smoothed at T>0, takes 
place when some parameters are varied and the pertinent stability condition is violated. 
In the present work we consider the model of a homogeneous Fermi liquid, which also 
displays a rearrangement of the ground state of the quasiparticle system, and investigate 
the scenario of the initial stage of the rearrangement. 
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2. Theoretical and computational aspects 



2.1. Two -connected Fermi sphere and the fermion condensate 

Let us start by recalling the Landau relation between the quasiparticle momentum dis- 
tribution n p (T) and the quasiparticle spectrum e p (T), 



E-, 



-i 



n p (r) = |l+exp ^' r ^' j , (1) 

(/i(T) is the chemical potential), which results from a variational equation 5F/5n p = 
(F is the free energy of the system) with th usual expression for the entropy [[17 



On one hand, eq. (1) is simply the Fermi-Dirac quasiparticle distribution over energies. 
On the other hand, this relationship is an equation for the quasiparticle distribution 
over the phase space. In fact, the quasiparticle energy is the variational derivative of 
the ground state energy functional E [n] with respect to the quasiparticle distribution, 
s p (T) = 5E /5n p (T), and, therefore, it is itself a functional of n p {T). 

It is postulated in the Landau theory that in a homogeneous isotropic Fermi liquid, 
like in a Fermi gas, the quasiparticle momentum distribution at T=0 has the shape of a 
fully occupied Fermi sphere rip\ P )=6(pF—p) (the maximum momentum pp is related to 
the density p of the system by the well known relation p=pp / '(37r 2 )) . The low-temperature 
behaviour of the quasiparticle spectrum corresponding to such a momentum distribution 



has the form [17]: 



e P (T)-^T)=i{p) + 0(T>). (2) 

The function increases monotonically in the vicinity of the Fermi momentum and 
changes its sign at p=pf- Its slope at this point, which is the group velocity of the 
quasiparticles on the Fermi surface vp=d^(p) / dp\ p=PF , is determined by one of the phe- 
nomeno logical parameters of the Fermi liquid theory, the effective mass M*= P p/vp. 

In a strongly correlated Fermi system, a quasiparticle momentum distribution min- 
imizing the energy functional £7o[n(p)] at T=0 may be located out of the corner point 
rip (p) of the functional space [n] and the low-temperature behaviour of the correspond- 
ing quasiparticle spectrum may differ from eq. (2). For instance, it was found in ref 



that quasiparticle energies which are equal to the chemical potential in a finite region of 
the momentum space can exist. In refs. ||, fLCfl , some model functionals E [n(p)] were 
introduced which, for certain values of the parameters, reach their absolute minimum for 
a momentum distribution characterized by a two-connected Fermi sphere 

np\p) = 9{pi-p) - 9(p 2 -p) + 0(p 3 -p). (3) 

A quite different quasiparticle ground state corresponds to systems with fermion con- 
densate [11-14]. Let us elucidate the main idea of the concept of a fermion condensate. A 
homogeneous and isotropic quasiparticle system with the fermion condensate is described 
by a singular solution of eq. (1) which corresponds to a quasiparticle spectrum linear in 
the temperature T within a finite region of momenta |TT|, |P2] : 

e p (T) - /i(T) = 7V (p) + o(T), Pi < P < Pf . (4) 
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Contrary to the Fermi liquid formula (2), there is no T- independent term in eq. (4). 
This means that at T=0 the quasiparticle spectrum has a plateau e p = \i in the region 
Pi<p<Pf. At T>0 the slope of the plateau is linear in T, and its position with respect 
to the chemical potential /x(T) is determined by the function vq(p) which is connected 
with the momentum distribution of quasiparticles in the condesate. Indeed, the singular 
solution of eq. (1), which can be easily obtained upon substitution the formula (4) into 
eq. (1), has the form n p (T)=n (p)+O (T) , where 

n {p) = {l+exp(i/ (p))} \ Pi<p<Pj (5) 

is the momentum distribution of the condensate quasiparticles at T=0. Outside the con- 
densate region, n (p)=l at p<Pi, and n (p)=0 at p>Pf |TT], 0. One can find the explicit 
form of n p (T) and e p (T) provided one knows the functional dependence e p (T)[n p (T)]. A 
set of model functionals was suggested in refs. [12-14, 16], each of them possessing the 
minimum at the singular solution, within well defined regions of functional parameters. 
In the present work we indicate the possibility of a scenario, in which the transition 
to the fermion condensate occurs through an intermediate structure, corresponding to a 
multiconnected quasiparticle momentum distribution. 

2.2. The effective functional 

The initial stage of the rearrangement, which is studied in this paper, is characterized 
by variations of the quasiparticle momentum distribution within a relatively thin layer 
around p^Pf (see results below). Under these conditions one can use the concept of 
effective functional, which is widely used in many-body theory. We consider the simple 
effective quasiparticle functional for the ground state energy: 

E [n p (T)} = J ^n p (T)dr+^Jv(p-p')n p (T)n pl (T)dTdr' (6) 
with the effective repulsive interaction 

V(P-P) = 7- J° , 3 - (7) 



The symbol dr in eq. (6) means integration over d 3 p'/(2n) 3 and summation over spin in- 
dices. Calculating the variatonal derivative of E with respect to quasiparticle distribution 
n p {T) one obtains the quasiparticle spectrum 

2 

e P (T) = {m + J y(p_p/) n ^ T) dT - (8) 

The functional dependence given by eqs.(7), (8) together with eq. (1) and the normaliza- 
tion condition 

n p (T)dr = p (9) 



are the set of equations to be solved for the quasiparticle disrtibution n p (T) and the 
spectrum e p (T). The value a=0.07pp was used in the calculations of the present work and 
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the behaviour of the system as the parameter Vq is varied was studied. The dimensionless 
parameter ■j=MVq/ (4:1t 2 pf) will be used in the discussion below. 

2.2. Numerical aspects 

Eq. (8) together with the expression (1) represents the nonlinear integral equation for 
the function s p {T). This equation was solved numerically by an iterative procedure with 
weighting factors. The 5-point Newton-Cotes quadrature formula with 5-point filter on 
output was used for numerical folding the distribution n p (T) with the effective interaction 
V(p,p'). The momentum grid had a step h p =5-10~ 5 pp. The accuracy of the numerical 
solution was determined upon substitution it into the initial equation. The permissible 
error, that is the maximum discrepancy between the left and the right hand sides of 
eq. (8), was fixed at 10~ 8 ep. The number of iterations necessary for reaching this accuracy 
is about 3 • 10 4 provided the iteration weight u>=0.001 is taken, which is optimal for a 
stability of the iterative procedure. It is worth to note that the results are independent 
of the point in the functional space taken as an initial one for the the iterative procedure. 
For example, the same solution for 7=0.50 at T=10~ 7 was obtained by starting from i) 
the solution for 7=0.50 at T=10~ 5 , ii) the solution for 7=0.48 at T=10~ 7 (here and below 
the temperature T is taken in units of e° F = p 2 F /2M). 



3. Results 

3.1. Instability of the Fermi-step momentum distribution 

We begin from estimating the value 7^°- ) of the interaction strength at which the necessary 
condition for stability of the system with the quasiparticle distribution rip (p) at T=0 is 



violated. This condition is fulfilled [11, O] provided the variation of the ground state 



energy Eq under any admissible variations of the distribution n(p), 

5E Q = J[ e (p)-v]5n(p)j^L, (10) 

is positive. Admissible variations of rip (p) are of the same sign as the difference p—pF, 
as dictated by Pauli principle. Hence, upon substitution of the energy e{pp) for the 
chemical potential /i in eq. (10) one can reformulate the necessary stability condition as 
a requirement for the value 

s(p)=2M £Mz£(|£) 

V ~Pp 

to be positive for each value of the momentum p |Tl|, |j"2f . Therefore, if the function s(p) 
has the first zero close to pp, the violation of the stability condition means that a bend 
appears in the curve e(p) in the vicinity of the Fermi momentum. For the distribution 
rip (p) the derivative de/dp can be easily calculated in the dimensionless form 

^(Q)(p) _ M de _ p | 7p F -f(p 2 +p 2 F +a 2 ) ^ {p+p F ) 2 +ct 2 _ 
Pf dp pp p Ap 2 (p—p F ) 2 +a 2 

In fig. 1, where the curves C (p) are shown for different values 7, one can see that the 
contact with zero takes place at the point p = p c ~ 0.97pF for 7=7^~0.415. Note that 
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the proximity of p c to pp justifies the substitution of the function s(p) with the function 
C^(p)- Therefore, the ground state with the quasiparticle momentum distribution rip (p) 
becomes unstable at 7>7^ ^ and its rearrangement takes place. 

3.2. Multiconnected momentum distribution 

To see how the ground state is arranged right beyond the transition point, let us look at the 
fig. 2, where the results of the calculations of n p (T) for 7=0.45 at different T are shown. At 
T=2-1CT 3 , the quasiparticle momentum distribution n p (T) has the shape characteristic of 
a system with fermion condensate (we shall discuss that in detail below). At T<2-1CT 3 a 
downfall appears in the distribution, which deepens as temperature decreases, and finally, 
at T=10~ 7 , one can hardly distinguish n p (T) from the two-connected Fermi sphere defined 
in eq. (3). The quasiparticle spectrum e p (T) corresponding to such distribution calculated 
at T=10~ 7 is shown in fig. 3. The spectrum of a fermion condensate has a shape with 
a plateau at a value exactly equal to the chemical potential \x at T=0 and getting little 
sloping at T>0 jllj |i~2| . Here we have a quite different behaviour, and the quasiparticle 
spectrum of the two-connected Fermi-like distribution equals \x at three points, which are 
the border pi of the inner sphere and the borders pi and p% of the spherical layer. The 
deviation of e p (T) from fi reaches the value ~ 2 • lO -4 ^ at the point of the minimum 
of the spectrum and the value ~ 2 • 10~ 6 ep at the point of its maximum. Despite the 
latter value is small, it is still higher than the estimated numerical error by two orders of 
magnitude. 

The behaviour of the two-connected Fermi sphere with increasing interaction strength 
7 is displayed in fig. 4. The spherical layer appearing beyond the transition point has, 
from the start, a finite thickness, while the gap between the layer and the inner sphere 
develops starting from a vanishing small width. The outer layer gets thicker and moves 
away from the inner sphere as the parameter 7 increases. What then happens with such 
a distribution? To understand that, let us investigate the stability of such Fermi sphere 
divided into two layers. We calculate the function C (p) 101 the momentum distribution 
in the form of eq. (3) and see when and where the critical change of the sign of that 
function occurs. An elementary calculation gives 

AD (ri ^ Md£ P , rv-ij^ 7(p 2 +P?+<* 2 ) , (p+ P i) 2 +a 2 \ 

C [P) ~Y F Tp-p~ F + U ' IT V (p-p^+ail {U} 

The function C (p) calculated for different values of 7 is displayed in fig. 4. The points 
of the maximum absolute values of the derivative dn/dp at T=10~ 7 were taken as the 
boundary momenta Pi,P2,P3 for that calculation. One can easily realize that inside the 
region 7^ < 7 < 7W ~ 0.452, the two points where C^(p) changes sign are located in 
such a way that the corresponding local minimum and maximum of e(p) lie in the domains 
where n(p) equals and 1 respectively. This means that the sign of the difference e(p)—fi 
coincides with that of the possible variations 5n(p) allowed by Pauli principle. Hence the 
considered distribution satisfies the stability condition. 

However one can see in fig. 4 that at 7 > 7W ~ 0.452, the situation changes. With the 
displayed behaviour of the function C^(p), there are regions where e(p)—fi>0, but with 
n(p) = 1. This means that the necessary stability condition is violated, since allowed 
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variations 5n(p) exist which lower the ground state energy. This results in the second 
rearrangement of the ground state of the system emerging at 7=7^ 1 - ) . It is shown in 
fig. 5 how the quasiparticle distribution is arranged beyond the second transition point 
7W. For definitness, the results of the calculations for 7=0.46 at different T are diplayed. 
The calculations show that the new layer appears for the above mentioned value of the 
coupling constant and the quasiparticle distribution n p (T) at T=10~ 7 is very close to the 
three-connected Fermi sphere 

(2) 

ny[p) = 9{pi-p) - 9{p 2 -p) + 0(p 3 -p) - 9(p A -p) + 6{p 5 -p). (14) 

The calculation shows that the appearence of new layers at increasing values of the pa- 
rameter 7 does not stop at the level of three-connected Fermi sphere. Figs. 6 and 7, where 
the quasiparticle distributions n p (T) are calculated for 7 = 0.48 and 0.50, show further 
divisions into layers of the momentum space. In particular, the distribution n p (T) for 
7=0.50 approaches a multi-connected Fermi sphere with lowering T, which is arranged as 
follows: the inner small occupied sphere with the radius ~ 0.85pF is surrounded by four 
spherical occupied layers with thickness ~ 0.3 — 0.4p^ divided by spherical empty layers 
with thickness ~ 0.1 — 0.2p F . The quasiparticle spectra e p (T) for 7=0.50 are shown in 
fig. 8. At T=10~ 7 , the spectrum crosses the line e=fi nine times, at the boundary of the 
inner sphere and at the boundaries of the spherical layers. 

4. Discussion 

4-1. Two scenarios of the rearrangement 

Thus beyond the first transition point 7^, the scenario of the rearrangement of the 
quasiparticle ground state at low T<10 -3 is the succession of transitions with increasing 
7. Each one of them results in the appearance of a new spherical layer in momentum 
space. Let us compare this scenario of the rearrangement and some distinctive features 
of the ground state under consideration with those of the fermion condensation. 

Let us recall the main features of the fermion condensation phenomenon. One of 
them is the plateau in the quasiparticle spectrum e p (T), with a value exactly equal to 
the chemical potential \x at T=0, in accordance with eq. (4), and which gets a finite slope 
at increasing temperature. Unlike systems with the fermion condensate, the spectrum 
e p (T=0) for the multi-connected Fermi-like distribution equals the chemical potential in 
a finite number of points, which are the boundaries of the spherical layers. The low- 
temperature expansion of such a spectrum has the form of eq. (2), typical of the Landau 
theory, with the non-monotone function changing its sign several times, unlike the 
monotone one for the usual Fermi liquid. At T=0, the ground state with the multi- 
connected Fermi-like quasiparticle distribution is not macroscopically degenerated, unlike 
the ground state of a fermion condensate. At the same time, there exist singularities of 
the density of states connected with the maxima and the minima of the function e p (T=0). 
These singularities gradually disappears with increasing T up to T ~ 2-10~ 3 , where the 
last twist of the spectrum corresponding to the outer layer is smoothed. At T>T the 
difference e p (T) —fi(T) becomes linear in temperature, like that for systems with a fermion 
condensate. 
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The other feature of systems with fermion condensate is the shape of the distribution 
n(p) given by eq. (5). Within the region occupied by the fermion condensate 0<n(p)<l, 
this corresponds to non-zero entropy of the fermion condensate at T=0. The contradiction 
with the Nernst theorem disappears providing correlations (e.g. superfluid) are taken 
into account, which immediately rearrange the ground state due to its degeneracy and 
re-establishes a zero entropy at T=0. The entropy of the state with the multi-connected 
Fermi-like quasiparticle momentum distribution equals zero at zero temperature because 
n(p) takes only the values and 1. The multi-layered distribution changes quickly with 
increasing T: the sharp boundaries of the layers are smoothed and the layers combine 
all together at T=Tq into a monotonically decreasing curve which is very similar to the 
momentum distribution of a system with the fermion condensate at that temperature. 
Being equal zero at T=0, the entropy of the system with the multi-connected Fermi-like 
distribution increases sharply with temperature due to quick smoothing of the function 
n p (T). The calculation show that at T=T the entropy reaches the value Sq^Qq/Q, the 
ratio of the phase volume of multi-layered region Q to that of the whole system Q. Just 
this value of the entropy would be characteristic at T~0 for a system with the fermion 
condensate occupying the phase volume Qq. At T>Tq the entropy becomes linear in 
temperature, like that of a system with a fermion condensate |12|, [16| . 

All these features of the momentum distribution, entropy, quasiparticle spectrum and 
density of states for a system with a multi-connected Fermi-like quasiparticle distribution, 
as well as the problem of validity of such a quasiparticle pattern, will be studied in detail 
in a separate paper. 

The scenario of the fermion condensation is characterized by the single critical value 7 C 
of the coupling constant, at which the fermion condensate arises in the system. The phase 
volume of the fermion condensate increases as 7 further increases, however the shape of 
the momentum distribution does not modify and no further qualitative changes occur 



Li], |12[ . The scenario of the rearrangement found for the model under consideration is 
different for different temperatures. At T=0 it is characterized by a sequence of critical 
values yj?\ each one corresponding to the appearence of a new ground state with a larger 
number of connectivity. The number of critical constants decreases as T increases, so that 
only one of them, 7^, survives at T>T . This means that for the considered model, the 
scenario of formation of the multi-connected Fermi sphere at T=0 gradually transforms 
into that of the fermion condensation with increasing temperature. 

4-2. Interplay between the multiconnected distribution and the fermion condensation 

Unfortunately, the computation time sharply increases with increasing the phase volume 
Qq of the layered distribution. This is the reason why the calculations of the present 
work are carried out up to the value of the interaction strength 7<0.5, corresponding 
to the initial stage of the rearrangement. What happens in the system at larger values 
of 7? To have a feeling of that, let us try to use the mechanical analogy, treating the 



momentum p as a spatial coordinate r \TQ [12|]. Then the problem of minimizing the ground 
state energy functional (6) at T=0 can be interpreted in the language of mechanics as 
a search for the statical equilibrium of the spatial distribution z/(r) of particles, which 
move in the external harmonic field U(r) = kr 2 /2 with the stiffness k=l/M and interact 
among each others by the repulsive force (7), the particle number being fixed by the 
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normalization condition (9). Unless the solution z/(r) of the mechanical problem exceeds 
2/ (27r) 3 , at least in one point, it can not be accepted as a solution n{p) because the latter 
should satisfy the additional restriction due to the Pauli principle. With strengthening 
interparticle repulsion, the considered mechanical system, obviously, expands and rarefies. 
As soon as the distribution u(r) is everywhere less than 2/(27r) 3 , one can conclude that 
it corresponds to the solution n{p) of the initial quantum problem. Since one can expect 
that the distribution v{r) should be smoothed and monotone, the fermion condensate 
seems to appear in the system at T=0 for large values of the interaction strength 7. This 
transition will be investigated in further publications. 

5. Conclusion 

In summary, the structure of the ground state of the homogeneous Fermi liquid with the 
strong interparticle repulsion is studied within the framework of the effective functional 
approach. The numerical investigation of the simple functional with strong repilsive effec- 
tive interaction showed that at fixed value of the interaction radius there exist the critical 
value of the interaction strength 7^, beyond which the ground state with the Fermi- 
step quasiparticle momentum distribution becomes unstable and the rearrangement of 
the ground state takes place. The scenario of the initial stage of that rearrangement with 
increasing 7 was found to be different for different regions of temperature. At T=0, there 
exist the set of critical constants 7W corresponding to the succession of transitions, each 
resulting in emerging the new spherical layer of the quasiparticle momentum distribution 
n(p)). The quasiparticle spectrum e(p) corresponding to such the layered distribution 
does not have the plateau, unlike the fermion condensate, and equals the chemical poten- 
tial in the finite number of points, which are the boundaries of the layers. The ground 
state with the multi-connected Fermi-like distribution possesses no macroscopic degenera- 
cy and the entropy of this state is zero at zero temperature. With increasing temperature 
the layers are quickly smoothed, so that at T~To~2-10~ 3 , there is no more reminiscence 
of the critical constants 7W with the exception of 7^. At T>T the scenario of the rear- 
rangement is that of the fermion condensation. The qualitative analysis showed that the 
found stucture with the multiconnected quasiparticle momentum distribution is the inter- 
mediate one, yielding the place for the fermion condensate with incresing the interaction 
strength. 
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Figure captions 



1. The function C (0) (p) calculated for 7 = 0.410, 0.415, 0.420. 

2. The quasiparticle momentum distributions n p (T) calculated for 7=0.45 at different 
temperatures. 

3. The quasiparticle spectrum (e p (T)—fj,)/e F calculated for 7=0.45 at T=10~ 7 . 

4. The quasiparticle momentum distributions n(p) and the function C^Kp) calculated for 
different values of the parameter 7. 

5. The same as for Fig. 2, for 7=0.46. 

6. The same as for Fig. 2, for 7=0.48. 

7. The same as for Fig. 2, for 7=0.50. 

8. The quasiparticle spectra {e p (T)—p)/e° F calculated for 7=0.50 at T=10~ 4 (solid line), 
T=10~ 5 (long dashes), T=1Q- 6 (short dashes), and T=10~ 7 (dots). 
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